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Abstract 



We study the effect of electron-electron interactions on the resonant- 
tunneling spectroscopy of the localized states in a barrier. Using a simple 
model of three localized states, we show that, due to the Coulomb interac- 
tions, a single state can give rise to two resonant peaks in the conductance 
as a function of gate voltage, G(V g ). We also demonstrate that an addi- 
tional higher-order resonance with V^-position in between these two peaks 
becomes possibile when interactions are taken into account. The correspond- 
ing resonant-tunneling process involves two-electron transitions. We have 
observed both these effects in GaAs transistor microstructures by studying 
the time evolution of three adjacent G(V g ) peaks caused by fluctuating occu- 
pation of an isolated impurity (modulator). The heights of the two stronger 
peaks exibit in-phase fluctuations. The phase of fluctuations of the smaller 
middle peak is opposite. The two stronger peaks have their origin in the same 
localized state, and the third one corresponds to a co-tunneling process. 

PACS Numbers: 73.20.Dx, 73.40.Gk, 73.40.Sx 
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Resonant tunneling through localized states in an insulating barrier (RT) has been ob- 
served in several types of semiconductor structures. It occurs when the energy of tunneling 
electrons is close to the energy of a localized state (LS). If a LS resides near the midpoint 
of the barrier its contribution to the conductance is of the order of e 2 /h. In the first ex- 
periments on diode structured the barrier represented a planar oxide layer sandwiched 
between two metallic contacts. In the later experimentsiHl a GaAs quantum well with LS's 
in it was separated from the contacts by two AlGaAs barriers. On small area barriers the 
authors were able to resolve the tunneling through individual LS's occuring at resonant bias 
voltages. 

Another group of experiments was carried out using lateral semiconductor structures^! 
with a two-dimensional electron gas. The barrier between the contacts is created by depleting 
electrostatically the region under the gate. The advantage of such a geometry is that the 
energy position of a LS with respect to the Fermi level in the contacts can be varied by 
changing the gate voltage V g . As V g is swept, different LS's in the barrier align with the Fermi 
level. As a result, narrow peaks are seen in the V g dependence of the ohmic conductance 
G. The fact that the shape of these peaks remains unchanged with decreasing temperature 
allows resonant tunneling to be distinguished from the phonon-assisted transportB. The 
magnitude of a conductance peak is determined by the well-known formula 

_e 2 4rTV 

^ max - h (ri + r r)2 > w 

where Ti and r r are the tunneling widths of the LS associated with electron escape into 
the left and right contacts respectively. Since T/, T r depend exponentially on the location 
of the LS within the barrier, the measurement of the G(V g ) dependence can be viewed as 
a resonant-tunneling spectroscopy of localized states in the barrier: the V g positions of the 
peaks reflect the energies of LS's, while their magnitudes characterize the spacial positions 
of LS's with respect to the center of a barrier. 

Such a spectroscopy implies that LS's are well separated in the direction along the 
barrier (y-direction in Fig. 1). Otherwise, the Coulomb interaction between different LS's 
becomes important. Due to this interaction, the energy position of a LS, as revealed in G(V g ) 
dependence, becomes a function of the occupation of the neighboring LS's. To illustrate this, 
consider two LS's with energies E\ and £2 (Fig. 1). Since E\ < £2 the level S\ is the first to 
"resonate" as the Fermi level in the contacts, /1, moves up with increasing V g . For /i > E\ the 
level Si becomes occupied. This causes an upward Coulomb shift, Z7i 2 , of the second LS. As 
a result the occupation number of the second LS will change from to 1 not at fi — £2 but 
at /j, = e 2 + Ui2- However, as long as the two-LS's configuration is considered, one can still 
attribute each peak in G(V g ) to an individual LS. In the present paper we demonstrate that, 
in general, this is not the case. Namely we will show that due to the Coulomb interactions a 
single LS can cause two peaks in G(V g ) dependence. We call this effect (resonant tunneling 
through the same LS at two different positions of /x) re-entrant resonant tunneling. We will 
then present the experimental evidence for this effect in GaAs transistor microstructures. 

Consider a configuration consisting of three LS's with energies £1,62 and £3 (Fig. 1(a)). 
Let us trace the evolution of their occupation numbers with /1 moving upwards. At fi — 
£1 = /ii the first LS becomes occupied. As already discussed, the second LS changes its 
occupation at \i = £ 2 + U\ 2 = Hi- Both charged LS's produce the Coulomb shifts, U23 
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and Uis of the third LS so that with further increasing /i the third LS gets charged at 
(j, = e 3 + U13 + U 23 = /i3- If, as it is shown in Fig. 1(a), only the first LS is close to the center 
of the barrier, while the other two are displaced towards the left and the right boundaries, 
one should expect one strong peak and two weak peaks in G(V g ). Note, however, that for 
above scenario to apply one should verify that after the charging of the third LS the energies 
of the first and the second LS still remain below the Fermi level. Indeed for p, > [13 these 
energies acquire positive shifts U\s and U23 respectively. For the first LS the corresponding 
condition yUi + t/12 + ^13 < ^3 can be rewritten as 

£ 3 - £1 > U 12 - U 23 . (2) 

This condition can be violated if the first and the third LS are close enough in energy. Then 
the question arises: what would the evolution of the occupation numbers be in this case? 
Note first that the meaning of inequality (|2|) is that the total energy of two electrons in the 
system of three LS's is minimal when they occupy the first and the second LS. This energy 
is equal to E\ + e 2 + U\ 2 - Consequently, the reversed inequality means that the minimum 
energy for the two electrons corresponds to the occupation of the second and the third LS 
and is equal to E 2 = e 2 + £3 + U 23 . If this is the case then the process of the charging of the 
system will go as follows. At [A — fii the first LS gets charged. At \i = £2 + £3 + £^23 — £1 = f* 
we have E\ + = E 2 . At this position of the Fermi level an electron from one of the contacts 
and an electron from the first LS occupy the second and the third LS, leaving the first LS 
empty. (Note that this position of the Fermi level is lower than that required to charge just 
the second LS, since /1* < fi 2 if the condition (|2|) is violated.) Finally, the first LS will be 
charged again at \i — S\ + U\ 2 + U\ 3 > ji* . 

From the evolution of the occupation numbers described above one should anticipate 
two strong resonant-tunneling peaks in G(V g ): at [A — E\ and at [A — E\ + U\ 2 + U 23 
corresponding to the resonance of the first LS. Note however, one more resonant-tunneling 
process is possible at an intermediate value of /x, namely at [A — //*. The "resonating" states 
are the states of two electrons. They are: 

(i) initial: an electron at the Fermi level in the left contact and electron on the first LS, 

(ii) intermediate: one electron on the second and another on the third LS, 

(iii) final: an electron at the Fermi level in the right contact and electron on the first LS. 
All three energies are equal at /i = //*. It can be shown that for such a process the magnitude 
of the peak in the conductance is still given by Eq. ([I]), but the values of T t and T r are now 
determined by two-electron transitions in the following way. 

One can write the general expressions for both widths: T; = 7r|VJ| 2 g, r r = 7r|K| 2 g, where 
g is the density of states in the contacts and Vi, V r are the corresponding matrix elements. 
In the case of the resonant tunneling through the first LS one has VJ = Voexp(— ■^ i ), V r = 
Voexp(— — ), where a is the localization radius, Ru, R\ r are the distances from the first LS 
to the left and to the right contacts respectively, and Vq is the prefactor. Then the formula 
(0) for G m ax takes the form 

Gmax = ^ COsh" 2 ( — - — \ . (3) 

For the two-electron transition between the states (i) and (ii) the matrix element V\ can 
be presented as 
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where Ri 2 is the distance from the second LS to the left contact and R\ 3 is the distance 
between the first and the third LS's. Two terms in the matrix element (|j) reflect the fact 
that two steps of the transition can be carried out in a different order (and thus through 
different virtual states). We emphasize that this matrix element is non-zero only due to the 
interactions. Indeed, it is easy to see that in the absence of the Coulomb shifts two terms 
in (|4]) cancel each other out. 

The matrix element V r of the transition between the states (ii) and (iii) has a form similar 

to m 



V r = Voexp 
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where R r3 is the distance from the third LS to the right contact and R\ 2 is the distance 
between the first and the second LS's. Using Eqs. (Q), (^j one can find the two-electron 
widths Ti = ir\Vi\ 2 g, T r = 7i\V r \ 2 g. Substituting Ti, T r into Eq. ([!]) we obtain the following 
expression for the maximal conductance 



G r 



— cosh 2 
h 



R12 + R13 — R12 — R 
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+ 7 



where the dimensionless constant 7 is defined as 

(C/12 - U 2 z){e 2 - ei)(e 2 - £1 + U 23 
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We see that the condition G max ~ 



and it can be 
is small. 



e 2 /h has the form R i2 + R 13 = R 12 + R r3 , 
satisfied even if the first LS is displaced from the center of the barrier, so that G 

The conduction involving virtual states considered here is similar to the elastic co- 
tunneling process in the Coulomb blockade^. In co-tunneling, electrons traverse a conduct- 
ing grain using virtual states when the transition via a real state is forbidden by the charging 
energy00. The main difference is in the number of virtual states. The quasi-continuum of 
states available in a metallic grain leads to the interference suppression of the co-tunneling 
rate0, while such a suppression obviously does not occur in the problem considered here. 

To study the effect of the re-entrant RT experimentally, we used sub-micron GaAs MES- 
FETs with lateral barriers00. Earlier, on these structures we demonstrated the presence of 
elementary RT channels in G(V g ) and I(V) characteristics and distinguished between one- 
and two-impurity RT. This lateral structure with a smooth barrier, Fig. 1, is well suited for 
testing the discussed model. It has an enhanced density of localized states near the contacts 
which will not give large peaks in G(V g ) unless two-impurity RT is involved^!, although a 
change in the charge of these states with varying V g can affect the states close to the middle 
of the barrier. By a standard lock-in technique we have measured the conductance of GaAs 
MESFETs with a donor concentration of 10 17 cm -3 . Highly doped source-drain contacts 
have metallic-like conduction. The barrier between the source and drain with localized 
states in it is produced by a gate of length L=0.2 /im in the direction of the current and 
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width W=20 /im in the perpendicular direction. The height of the barrier is controlled by 
the gate voltage which shifts the localized states in the barrier with respect to the Fermi 
level in the contacts. 

The experimental demonstration of the re-entrant resonant tunneling should provide an 
evidence that the two peaks in G(V g ) have their origin in the same localized state. Firstly, 
we expect that the heights of these peaks should be similar. (A possible difference in the 
heights can be due to some effect on the values of T; and r r in (1) when levels 2 and 3 are 
charged.) Secondly, the external fields should influence both peaks in the same way, if the 
peaks originate from the same energy level. If this field is a random electric field caused 
by the occupancy fluctuations of an isolated impurity (modulator), such a modulator can 
serve as a local probe to the arrangement of localized states in the barrier. To explore this 
possibility we have measured G(V g ) curves repeatedly for a long period of time (several 
hours) and analysed the fluctuations of peaks in G(V g ) with time. 

Fig. 2 shows an example of G(V g ) at T=100 niK with the peaks which we associate with 
one-impurity RT channels. Among a dozen peaks there is a group of three neighbouring 
peaks labelled A, B and C that have two distinct states at different times. Two-level con- 
ductance fluctuations in mesoscopic structures (Random Telegraph Signal) are commonly 
connected with a single impurity placed close to a conducting path of small size0. The fluc- 
tuations of the occupancy of this impurity produce a random electric field which modulates 
the conductance of the path. The surprising feature of the fluctuations of the three peaks 
in Fig. 2 is that these fluctuations are coherent: the decrease of peaks A and C and increase 
of peak B occur simultaneously. Fig. 3(a) shows the fluctuations of the peak amplitudes 
with time. Comparison of the relative magnitudes of these fluctuations shows that for the 
peaks A and C they are very close: Ghi g h/Gi ow is 1.52 for peak A and 1.58 for peak C. 
The fluctuations in the amplitude are accompanied by distinct but less pronounced shifts 
in the positions of the peaks, Fig. 3(b). We have also found a similar group of three peaks 
with coherent fluctuations on another sample from the same wafer. The only difference in 
the second sample is that the fluctuations of the middle peak B, which has the smallest 
amplitude, are hardly detectable. 

The fact that the observed fluctuations of the peak amplitude are simultaneous and have 
the same relative magnitude suggests that the peaks A and C can be considered as being 
due to the same LS (state 1 in our model). Let us examine first an alternative explanation 
where these peaks are due to different LSs. We will see that this straightforward view can 
be dismissed as highly improbable. Expression (1) for the tunneling peak amplitude can be 
reduced to G max = (e 2 /h)(T min /T max ) = (e 2 /h) exp[-(R max - R min )/a], where T min (T max ) 
is the minimal (maximal) value of Y h Y r and R ma x(Rmin) is the corresponding distance 
from the LS to the contacts, Fig.l,b. Similar values of G max for the two peaks means that 
the other LS has to lie on the same line along the barrier as state 1 (or on the symmetric 
about the centre line). The equal effect of the modulator on the two states means that their 
distances to the modulator have to be the same. Thus, for a given state 1 there is only one 
possible place for the other state to satisfy the experimental facts: at the intercept of the 
sphere with the modulator in its center, with the line in the 2D plane of the barrier, see 
Fig. 1(b). The uncertainty in the position of this state in the plane is of the order of a 2 , 
where a is the localization radius, and in the energy scale it is of the order of T max . With 
an accuracy of a numerical factor, the probability to meet such a localized state P]_ 2 is 
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g2D0 2 T maX) where g 2 D is the two-dimensional density of states. In our structures the value 
of g2D is about ten times less than the conduction band value 3 x 10 10 cm _2 meV _1 , and 
a ~ 100 A. The rate of the Fermi level shift with gate voltage a = dfi/deV g ~ O.2HH gives 
Fmax ~ 0.2 meV from the width of the peaks in the V g — scale. For the probability P1-2 we 
then obtain a very small value of 6 x 10~ 4 . 

For the re-entrant tunneling view we have to estimate the probability to meet two states 
in the vicinity of the first LS. These states should lie within an area of L c xL c for the Coulomb 
interaction not to be screened by the metallic gate, where L c is the distance between the 
gate and the conducting channel ( in our structures L c « 1000 A and is comparable to an 
effective barrier length ). The characteristic energy interval of these two states should be 
now of the order of the Coulomb energy Uj ~ e 2 /nL c rs 1 meV. This value agrees with the 
estimation of Uj from the separation AV g between peaks A and C in Fig. 2, which in the 
re-entrant tunneling model is Uu + U13 = aeAV g ~ 1 meV. The probability P1-23 to meet 
the two impurities is then (g2DUijL 2 c ) 2 « 0.1. The latter is two orders of magnitude larger 
than the value of P1-2 and it is large enough to make the effect of re-entrant tunneling 
observable. 

In our explanation the only feature of the modulator that is important is that it has a 
local effect on the surrounding localized states. Its most probable position is between the 
gate and the conducting channel as is indicated in Fig. 1(b). This conclusion is drawn from 
the fact that we do not see incoherence in the fluctuations of the peaks A and C which 
are shifted in V^-scale by AV g =6 mV. It means that this change in V g does not affect the 
average occupancy of the modulator. For an energy level under the gate such a change in V g 
will indeed have a small effect, as its occupancy is formed under the condition of a strong 
electric field between the gate and channel, V g « — 1.7 V. The fluctuations of the modulator 
charge in time can result from a small gate-channel leakage current, as it was seen inS 

We explain the difference in the amplitudes of peaks A and C by the influence of charged 
levels 2 and 3 on T min and T max in the expression for the first level resonance G max = 
(e 2 /h){T min /T max ). It is interesting to note that the modulator produces a similar change 
in the peak amplitudes but without significant change in the V g - position of the peaks (the 
shift accompanying the amplitude modulation is ten times less than the separation between 
peaks A and C). Within our model, with the second and the third LSs positioned on the 
opposite sides of the first LS, this can be understood as follows. Although states 2 and 3 
are closer to the state 1 than the modulator and thus produce a stronger Coulomb shift 
of the energy position of state 1, their effects on T min and T max are partially compensated 
in G max which is the ratio of the two values. The modulator, however, is most probably 
located assymetrically with respect to the first LS and it mainly affects one of the two Ts, 
so that the compensation does not occur. 

Peak B in Fig. 2 has a much smaller amplitude than peaks A and C. In addition, its 
modulation is in anti-phase with that of the peaks A and C: the decrease of the peak B 
corresponds to the increase of the other two. This implies that its origin is not in the state 1. 
For the states 2 and 3 we would normally expect to see two small peaks, unless these states 
are positioned too far from the middle of the barrier to give a measurable conductance. The 
presence of only one resonance allows us to suggest that the peak B in Fig. 2 is due to co- 
tunneling with all states 1, 2 and 3 involved. Note that for peaks A and C the modulation 
results from the change of the ratio Ti/T r in (1) by the electric field of the modulator, 
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whereas for co-tunneling (6) it is the ratio Ti/T r that determines the peak amplitude. It 
is quite plausible that when one ratio increases the other decreases and vice versa, which 
would lead to the anti-phase modulation of the co-tunneling peak. 

In conclusion, we have demonstrated that charging effects can drastically change the 
resonant-tunneling spectroscopy of localized states in the barrier. These effects cause the 
doubling of the single-impurity peaks and also give rise to additional peaks due to many- 
electron tunneling. 

The work at the University of Minnesota was sponsored by NSF Grant DMR-9423244. 
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FIGURES 



FIG. 1. a) Schematic diagram of the barrier with resonant tunneling through the first localized 
state (states 2 and 3 are not occupied), b) Spacial position of the localized states and a modulator 
in the potential barrier under the gate. The part with the state 1 and modulator M illustrates 
the estimation of the probability for the peak C to be due to an independent localized state. The 
part with three localized states (1,2 and 3) corresponds to the probability to meet three interacting 
states. 

FIG. 2. Two consecutive measurements of the ohmic conductance as a function of gate voltage. 
Peaks A,B and C show fluctuations between distinct two states while the neighbouring peaks are 
time independent. 

FIG. 3. a) Fluctuations of the amplitude of peaks A,B and C with time, b) Fluctuations of 
the positions of these peaks in V^-scale accompaning the amplitude fluctuations. The peak D at 
V g = —1.734 V is the neighbour of peak C and it does not show coherent fluctuations. 
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Fig.l 
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